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£> Abstract 

Absent any indirect tests on the thermal history of the Universe prior to the formation of light 
nuclear elements, it is legitimate to investigate situations where, before nucleosyntheis, the sound 
r^. , speed of the plasma was larger than c/y/3, at most equalling the speed of light c. In this plausible 

extension of the current cosmological paradigm, hereby dubbed Tensor-ACDM (i.e. TACDM) 

OO 

scenario, high-frequency gravitons are copiously produced. Without conflicting with the bounds on 

the tensor to scalar ratio stemming from the combined analysis of the three standard cosmological 

data sets (i.e. cosmic microwave background anisotropies, large-scale structure and supenovae), the 

spectral energy density of the relic gravitons in the TACDM scenario can be potentially observable 

by wide-band interferometers (in their advanced version) operating in a frequency window which 

ranges between few Hz and few kHz. 
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The only direct informations on the early thermal history of the Universe come, at present, from a 
background of relic photons which last scattered the electrons at an approximate redshift of z^ cc ~ 1087 
according to the 5-yr WMAP data release [1, 2]. The scrutiny of the Cosmic Microwave Background 
(CMB) observables is always conducted within a commonly accepted framework, i.e. the so-called 
ACDM paradigm where A qualifies the dark-energy component (parametrized in terms of a cosmolog- 
ical constant) and CDM qualifies the (cold) dark matter component. The ACDM scenario represents 
a useful compromise between the available data and the number of ascertainable parameters. A class 
of plausible completions of the ACDM model contemplates the addition of a post-inflationary phase 
expanding at a rate which is slower than radiation. From the point of view of the fluid properties, the 
sources generating such a dynamics are often called stiff. The spectral energy density of the gravitons 
reentering the Hubble radius during the stiff phase increases with the frequency rather than being 
nearly constant as in the conventional ACDM paradigm. Such an extension requires two parame- 
ters: a typical frequency scale, be it v s (corresponding to the end of the stiff epoch) and the slope of 
the spectral energy density during the stiff phase. The supplementary parameters characterizing this 
scenario (which will be dubbed, in what follows, as TACDM for tensor-ACDM) can be determined 
by analyzing the three conventional cosmological data sets (i.e. CMB [1, 2], large-scale structure 
[3, 4] and supernovae [5, 6]) in conjunction with the forthcoming data of wide-band intereferometers 
[7, 8, 9, 10]. At the moment interferometers are only able to provide interesting upper limits on the 
spectral energy density of the relic gravitons [11]. The foreseen sensitivities of the so-called advaced 
Ligo [7] will still be inadequate to probe the relic gravitons produced within the conventional ACDM 
scenario. Nonetheless the very same sensitivities of the interferometers in their advanced version will 
be definitely sufficient to probe directly the parameter space of the TACDM scenario. 

Consider therefore the evolution of the tensor modes in conformally flat background geometries 
which are, incidentally, the ones currently preferred in the context of the ACDM paradigm [1, 2]. A 
conformally flat background geometry in four space-time dimensions, by definition, is characterized 
by a metric ~g„ u = a 2 ^)^^ where rj^ is the Minkowski metric with signature mostly minus and r 
is the so-called conformal time coordinate. The tensor fluctuations of the geometry are defined with 
respect to the three-dimensional Eucledian sub-manifold as 

5^g l0 = -a 2 h l0 , ,5 t ( V = -^, 4 V = "^, dih) = h\ = 0, (1) 

where Latin indices run over the spatial dimensions; <5j: and 5\ denote, respectively, the first and 
second order tensor fluctuations of the corresponding quantity. Since hij is a (divergenceless and 
traceless) rank-two tensor in three spatial dimensions, it carries two physical polarizations. Defining 
three mutually orthogonal directions as k\ = ki/\k\, rhi = m.i/\fh\ and h{ = rij/|n|, the two polarizations 
of the gravitons in a conformally flat background are nothing but 

€ if (*0 = (mimj ~ fiifij), e-if (k) = (mhj + fUrhj), e ij 4j = 25 aa'- (2) 



By perturbing the Einstein-Hilbert action to second order in the tensor amplitude hij, the action for 
the gravitons can be written, up to total derivatives, as 

S gw = 5 ( ^S = ±j d 4 *,/^ r u dphijdvW, £ P = V8^G = ^- = ±-. (3) 
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Up to a rescaling of the amplitude in terms of the Planck length, the canonical normal modes of the 
action (3) are given by pjj = ahij since, in a conformally flat background, \/—g g^ u — ► a 2 (T)rj^ u . The 
mode expansion of the canonical field operator is thus given by: 
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where [or , , oL\i\ = 5' 3 ' {k—p)8\\> . It will be hereby assumed that the field operators are in the vacuum 
at the onset of the inflationary evolution. Thus the initial state |0) (annihilated by at ,) minimizes 
the tensor Hamiltonian when all the wavelengths of the field are shorter than the event horizon at the 
onset of the inflationary evolution (see, for instance, [12]). In Eq. (4) f^\ are the (complex) tensor 
mode functions obeying 

fk,x = 9k,x, g'k,x = -[k 2 -{n' + n 2 )]h,x, n = ^, (5) 

where the prime denotes a derivation with respect to the conformal time coordinate r. Defining p t 
and p t as the total pressure and as the total energy density of the plasma, the Friedmann-Lemaitre 
equations read: 

H 2 = ^-a 2 p u n 2 -n' = 47TGa 2 ( Pt + Pt ), p( + 3H(p t +Pt)=0. (6) 

If the inflationary phase is suddenly followed by the radiation-dominated phase, the energy-density of 
the inflaton is instantaneously converted into a radiation. This approximation is customarily employed 
to assess the number of inflationary e-folds [12, 13]. Given our ignorance on the thermal history of the 
plasma prior to nucleosynthesis, the inflationary phase might not be suddenly followed by the radiation 
dominated phase [14, 15]. Provided the transition between inflation and radiation is sufficiently stiff 
(and long) high-frequency gravitons can be copiously produced [14, 15]. A relativistic plasma is said 
to be stiff if its sound speed is larger than the sound speed of a gas of ultra-relativistic particles 2 i.e. 
l/\/3. The total sound speed and the barotropic index are defined, respectively, as: 
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where, in the second equality defining c 2 t , Eq. (6) has been used. In the primeval plasma, stiff phases 
can arise: this idea goes back to the pioneering suggestions of Zeldovich [16] in connection with the 



2 Natural units h = c = &b = 1 are used throughout the script. 



entropy problem. If an inflationary phase precedes a stiff phase the spectral energy density of the 
relic gravitons increases with frequency and the typical length of the stiff epoch can be determined 
by back-reaction effects [14]. In [15] the techniques of [14] were applied to assess the spectral energy 
density in the models of quintessential inflation which were developed in [17]. There were various 
reprises of these ideas (see, for instance, [18] and references therein). A (causal) upper limit on wt and 
c st is the speed of light, i.e. Wt < c st < 1 [19]. 

Collisionless species couple to the tensor modes of the geometry. Defining as Ily the anisotropic 
stress of the plasma we will actually have that below temperatures O(MeV), i.e. after weak interactions 
fall out of thermal equilibrium, the evolution equations for the classical amplitude corresponding to 
the quantum operators of Eq. (4) reads 

Hij - V 2 Liij - (H f + n 2 )fHj = -l6irGa 3 U i:j . (8) 

The coupling to the anisotropic stress induces computable differences on the spectral energy density of 
the relic gravitons. The effects of neutrino free streaming has been investigated both semi-analytically 
and numerically in [20] (see also [21, 22, 23]). With this caveat on collisionless species, Eqs. (5)-(6) 
can be solved numerically; the spectral energy density of the relic gravitons (and the related power 
spectrum) can then be assessed. 

The definition of the energy-momentum pseudo-tensor of the gravitational field always involves 
a certain degree of ambiguity. After getting rid of the tensor structure by making explicit the two 
physical polarizations, the action of Eq. (3) is just the action of two minimally coupled scalar fields 
in a conformally flat geometry of Friedmann-Robertson- Walker (FRW) type. The energy-momentum 
pseudo-tensor of relic gravitons in a FRW background just given by [24] 
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where the second equality follows from the first by using hij = J2x h(\) e ij an d by recalling the orthog- 
onality condition appearing in Eq. (2). In a complementary perspective [25], the energy-momentum 
pseudo-tensor is instead defined from the second-order fluctuations of the Einstein tensor, i.e. 

r ; = -^t (2) ^, Ql = r; - ^s;r, (10) 

where the superscript at the right hand side denotes the second-order fluctuation of the corresponding 
quantity while the subscript refers to the tensor nature of the fluctuations. The two definitions seem 
very different but the energy densities and pressures derived in the two approaches give coincident 
results as soon as the corresponding wavelengths are inside the Hubble radius, i.e. k > TC. In the 
opposite limit Eqs. (9)— (10) seem superficially different but give consistent quantitative results once 
they are compared on a particular background geometry [26]. 



By definition, pgw(x,t) = (0|Tq(x, t)|0) where |0) is, again, the state annihilated by dk t \- Re- 
calling that the mode functions of each polarization coincide ( i.e., in Eq. (5), fk,® = fk,® = fk an d 
analogously for g^) the spectral energy density in critical units can then be expressed as: 
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The spectral energy density of the relic gravitons can be related to the power spectrum which is, by 
definition, the Fourier transform of the two-point function evaluated at equal times, i.e. using Eq. (4) 



(0|^(x,r)^'(y,r)|0) = J d\nkP T (k,rf-^, V T (k,r) = 4£l-^-\f k (r)\ 2 , (12) 

where r = \x — y\. Quantum fluctuations present during the inflationary phase are amplified with 
nearly scale-invariant slope. The inflationary power spectra are then parametrized in terms of the 
tensor and scalar spectral indices, i.e., respectively, nx and n s : 

Arr / k \ "T / k \ ns_1 

„_£, V T (k)^A T (-) , V K(M= M(-) , (13) 

where k p = 0.002 Mpc~ is the so-called pivot wave-number which corresponds to an effective multipole 
£ e g ~ 30. In the context of the ACDM paradigm, the 5-yr WMAP data alone imply An = 2.41 x 10 -9 
(slightly different values can be obtained if different data sets are combined but these differences do 
not affect the features addressed here). The tensor amplitude is therefore estimated by setting limits 
on tt which is, by definition, the ratio between the tensor and the scalar amplitudes evaluated at the 
pivot scale k p . The inferred upper bounds on n, range from rx < 0.2 (in the case of the WMAP 5-yr 
data alone [1, 2]) up to rx < 0.43 when the WMAP 5-yr data are combined with the large-scale [3, 4] 
and supernova data [5, 6] (see also the thorough analyses reported in [1, 2]). In the minimal version 
of the inflationary dynamics the tensor spectral slope (i.e. nx), the slow roll parameter e as well as rx 
are all related: 

n T ~-y~-2e, e= __> , (14) 

where e measures, as indicated, the (slight) decrease of the Hubble H rate during the quasi-de Sitter 
phase of expansion and the overdot in the last equation denotes a derivation with respect to the cosmic 
time coordinate. Usually Eq. (12) is computed at the present epoch and then, in a second step, the 
spectral energy density of the relic gravitons is derived [27] (see also, for instance, [21, 22, 23]). The 
spectral energy density can be also directly assessed by numerical means without passing through the 
transfer function of the amplitude: this will be the approach followed here. Within the first strategy 
the power spectrum is given by 
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where 3 , according to [21], c\ = 1.34 and b\ = 2.50. In Eq. (15) ji(y) = (siny/y 2 — cos y/y) is the 
spherical Bessel function of first kind which is related to the approximate solution of the evolution 
equations for the tensor mode functions whenever the solutions are computed deep in the matter- 
dominated phase (i.e. a(r) ~ r 2 ). To obtain the spectral energy density, Eq. (11) must then be 
evaluated in the limit k 2 3> TL 2 (i.e. wavelengths inside the Hubble radius). In the latter limit the 
tensor mode functions satisfy |<?a;(t)| ~ kfk(r) and Eq. (11) then gives: 

k 2 36i / k \ ™ T 

^Mk,ro) = T ^V T (k,ro), fc ^ q ° Gw(fc ' T0) " 8a^ V ^ q l^J " (16) 

Since V"c{k,To) oscillates also ^gwC^to) wm oscillate. In the limit k 3> k cq the cosine will dominate 
the expression of ji(fcro)and the second result of Eq. (16) arises by replacingcos 2 (/cTo) — ► 1/2. If we 
take b\ = 2.5 in the second relation of Eq. (16), then 3&i/8 = 15/16 = 0.9375. If we take instead our 
results (i.e. b\ = 2.683) we will get, for the same quantity, 1.006. What appears in Eq. (15) is the 
transfer function of the tensor amplitude which literally transfers the power spectrum Vt inside the 
Hubble radius. 

In a complementary perspective, the consistent numerical solution of Eqs. (5)-(6) allows for a 
numerical calculation of Qc^/(k, r) according to Eq. (11). Instead of fitting the final result in terms of 
a putative (semi-analytic) amplitude for the mode function, the momentum (or frequency) profile of the 
spectral energy density will be obtained directly by numerical methods. As usual, initial conditions for 
the numerical integration are given for fcr<l. The system is then followed through Hubble crossing 
(i.e. kr ~ 1). Finally, when fer> 1 the expression of A p (/c,r) can be read-off in the asymptotic 
regime. In Fig. 1 (plot at the left) the numerical integration across the radiation-matter transition 
is illustrated. Instead of phrasing the numerical integration in terms of k and r, it is practical to use 
x = kr and k = k/k cq as preferred variables. To be accurate on the initial conditions a fully analytic 
solution of Eq. (6), valid across the radiation- maatter transition, can be safely employed: 
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In the limit x <C 1 the initial conditions for the mode functions are determined directly (and up to 

phase factors) from Eq. (13). Since the system is linear, the tensor mode functions can be always 

rescaled through their initial value; the energy transfer function is therefore defined by the following 

limit 

limApfoz) =T 2 (n)A p (K,Xi), x\ < 1. (18) 



3 By repeating the analysis of [21] we obtained cii = 1.260 and b\ — 2.683 which is fully compatible with the results 
of [21]. In the approach of [21] the calculation of the amplitude transfer function, in fact, involve a delicate matching on 
the phases of the tensor mode functions. Conversely, if the transfer function is computed directly for the spectral energy 
density, the oscillatory contributions are suppressed as the wavelengths get shorter than the Hubble radius (see below). 
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Figure 1: The spectral energy density (see Eq. (11)) is integrated across the radiation matter transition 
for different values of k = k/k cq (plot at the left). A finer grid in k (plot at the right), allows for the 
computation of the energy transfer function whose form can be fitted by an approptiate analytical 
expression (see Eq. (19). 



In Fig. 1 the results of the numerical integration are reported in terms of A p (n,x) for different values 
of k (see plot at the left). Always in Fig. 1 (plot at the right), Tjj (k) can be computed numerically: 
the diamonds correspond to the numerical points and the full line (in the plot at the right) is the 
numerical fit obtained by means of standard methods in the analysis of the regressions: 
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For a successful numerical determination of T p (n) the initial integration variable should be sufficiently 
small (i.e. x\ = kr\ <C 1) in such a way that, at the initial time, the mode k eq = r" 1 had a corresponding 
wavelength much smaller than the Hubble radius at t;. Second, Xf should be sufficiently large so that, 
effectively, A p (xf,n) is constant up to terms 0(l/x{) (see also below Eq. (22)). Finally, the grid in 
k should be sufficiently fine to allow for a reasonable fit. Using Eq. (19), the spectral energy density 
can be written, in the absence of free streaming, as 
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where (3 = 6.33 has been determined numerically assuming a smooth transition between inflation and 
radiation [28]. Equation (20), unlike Eqs. (15)-(16), is not strongly oscillating. The rationale for this 
difference is that, when computing A p (n, x), the oscillating contributions get dynamically suppressed 
as the wavelengths get shorter than the Hubble radius. A way of understanding this effect is to notice 
that the crudest approximation for the mode functions in the limit kr ^> 1 are simple plane waves, 
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Figure 2: The spectral energy density is integrated across the transition between the stiff epoch and 
the radiation-dominated epoch for different values of k (plot at the left). Following the same procedure 
used in the case of Fig. 1 the energy transfer function can be obtained and fitted by the appropriate 
analytic expression (plot at the right). 



i.e. 



fk(r) 



c + {k)e~ tkT + c-(k)e 



ll-.T 



9k(r) 



-°V2 



c + {k)e 



-ikr 



C-(k)e 



ikr 



Using Eq. (21) into Eq. (11) and enforcing the limit x — ► Xf S> 1, 

A p (K, Xi ) = k(\c + (k)\ 2 + |c_( K )| 2 ) + oQ-\ 



(21) 
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which proofs that the oscillating contributions are suppressed and that A p (k, Xf) is proportional to 
what are called, in the jargon, mixing coefficients. The considerations developed in the case of the 
radiation-matter transition also apply, for instance, to the stiff-radiation transition. In Fig. 2, for 
instance, the transition between a radiation-dominated phase and a stiff phase (with u> t = 1) is 
illustrated. This time the energy transfer function will be increasing with the wavenumber (see Fig. 
2, plot at the right) and the energy transfer function will be given, this time, by 

/ h \ V 4 / h \ V 2 / k\ / k\ 

T 2 (k/k s ) = 1.0 + 0.204 f — J - 0.980 f — J + 3.389f — J - 0.067 f — J In 2 (k/k s ), (23) 



where k B = t^T 1 and r s the time at which the plasma becomes dominated by radiation. The fact that 
the spectral energy density increases linearly (up to logarithmic corrections) fits with the analytical 
results of [14, 15] where, however, the slow-roll corrections were neglected. Further details on this 
approach will be given in a forthcoming paper [28]. 

The outlined computational procedure allows for a reasonably accurate estimate of the spectral 
energy density of the relic gravitons in a variety of models. In Figs. 3 and 4 the spectral energy density 
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Figure 3: The spectral energy density of the relic gravitons is illustrated in the case of the (conven- 
tional) ACDM paradigm supplemented by the tensor to scalar ratio rr- The parameters are fixed to 
the best-fit values derived by comparing the ACDM paradigm with the WMAP 5-yr alone [1, 2]. 



is reported, respectively, in the conventional case and in the TACDM scenario taking into account, in 
both cases, the late-time effects which can marginally reduce the amplitude. Depending upon R v (i.e. 
the neutrino fraction in the radiation plasma), the tensor amplitude and the spectral energy density 
get reduced. For three families of massless neutrinos (as implied by the WMAP 5-yr best fits and as 
assumed in the pivotal ACDM paradigm) R v = 0.405 and the amount of suppression is, approximately, 
0.64 of the value VIq^{v,tq) has when the very same effect is not taken into account. 

The effect of a progressive reduction of relativistic degrees of freedom has been approximately 
taken into account. In the least favourable case the reduction of the relativistic degrees of freedom is 
flat in frequency and proportional to (#p/<?po)(s , s/.S f so) 3 where <? p o = 3.36,<? s o = 3.90 [21, 22]. Note 
that g p and g s are the relativistic degrees of freedom appearing, respectively, in the energy and in the 
entropy density. Finally, there is a modification in the spectrum connected with the late dominance 
of the dark energy [21]. The most prominent effect is independent on the frequency: the spectral 
energy density is suppressed by an extra-factor, i.e. (Omo/^a) 2 - In the case of the WMAP 5-yr data 
alone, the ACDM paradigm gives Oa = 0.742 and f^MO = 0.258. Intuitively this means that that M p 
(appearing in Eq. (20)) is further suppressed by a factor 0(0.120). In Fig. 3 (plot at the left) /iq^gw 
is illustrated as a function of the frequency u = k/(2ir) by taking into account all the late-time effects 
mentioned above. The pivot frequency u p = 3.092 aHz corresponds 4 to the pivot wavenumber of Eq. 
(13). The spectral energy density (see Fig. 1 plot at the left) consists of a decreasing region (at low 



4 Whenever needed, the prefixes of the International System of units will be consistently adopted: laHz = 10 18 Hz, 
lfHz = lCT 15 Hz and so on. 



frequencies) which is followed by a nearly scale-invariant plateau for frequencies v > i/ eq where 

, eq = *s = i.28i x icr 17 f M^E) f ^"hd \ - 1/2 (24) 

q 2vr V0.13267 V4.15 x 10- 5 / v 7 

is the frequency corresponding to matter-radiation equality 5 . The WMAP 5-yr collaboration [1, 2] 
give an experimental determination of k cq , (i.e. k cq = 0.00999t 'ooo27 Mpc -1 ) which is fully compatible 
with the analytical estimate of Eq. (24). According to Fig. 3, /iq^gw(^) r o) decreases exponentially 
for v > f max where 

While f cq does not depend upon the specific model, u max depends, in principle, from the amount of 
redshift between the end of inflation and the present epoch. The shallow depression arising in the 
nearly scale-invariant plateau of Fig. 3 (plot at the left) is due to neutrino free streaming and it is 
present for z/ e q < v < z^bbn where 

- - -^ >< i »- 1 (^0 1/4 (M (ife) 1/4 »* - °- 01 -*■ < 2 « 

The frequency band of the terrestrial interferometers [8, 7, 9, 10] ranges between few Hz and 10 kHz 
with a maximum in the sensitivity to a stochastic background 6 for, approximately, v^M — 0.1kHz. 
Since f cq < u^m < v max , Fig. 1 implies (plot at the left) that /iq^gw (^lv , To) — 10~ 17 . To be 
even more quantitative, in Fig. 3 (plot at the right), /io^Gw( z/ LV ) r o) is illustrated as a function 
of JT- In the same plot, the dot-dashed curve refers to the standard case discussed in Eq. (14); 
the full and dashed curves refer instead to the situation where the spectral index depends upon the 
frequency as »x = — ^T/8 + (rT/16)[(n s — l) + (rx/8)] In (z//i/ p ). Figure 3 shows that, in both situations, 
^o^Gw( l/ LV)T"o) — OIO^ 17 given the current limits on rx- 

In the case of an exactly scale invariant spectrum the correlation of the two (coaligned) LIGO 
detectors with central corner stations in Livingston (Lousiana) and in Hanford (Washington) might 
reach a sensitivity to a flat spectrum which is [29] 

hi ^gw(^lv) - 6.5 x 10~ n {^Y-\ SNR 2 , z^ v = 0.1 kHz (27) 

where T denotes the observation time and SNR is the signal to noise ratio. Equation (27) is in close 
agreement with the sensitivity of the advanced Ligo apparatus [7] to an exactly scale-invariant spectral 



In Eq. (24) f^MO and £2ro are, respectively, the critical fractions of matter and radiation of the putative ACDM 
model. 

6 The sensitivity to a given signal depends upon various factors. For intermediate frequency the signal to noise ratio 
is also sensitive to the form of the overlap reduction function which depends upon the mutual position and relative 
orientations of the interferometers. The overlap reduction function effectively cuts-off the integral which defines the 
signal to noise ratio for a typical frequency v ~ l/(2d) where d is the separation between the two detectors. 
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Figure 4: The spectral energy density of the relic gravitons in the case of the TACDM scenario. The 
parameters are fixed to the best-fit values derived by comparing the ACDM paradigm with the WMAP 
5-yr alone [1, 2]. 

energy density [30]. Equation (27) together with the plots at the right in Fig. 3 suggest that the 
relic graviton background predicted by the ACDM paradigm is not directly observable by wide-band 
interferometers in their advanced incarnation. The minuteness of fo(j^Gw( i/ LV>' 7 o) stems directly from 
the assumption that the inflationary phase is suddenly followed by the radiation-dominated phase. 

Let us then posit that between the end of inflation and the onset of the radiation-dominated phase 
a sufficiently long stiff phase takes place. In this case the spectral energy density of the relic gravitons 
will increase for frequencies larger than u$ = k s /(2ir). Assuming that the inflationary phase is be of 
quasi-de Sitter type and characterized by a given value of rx, it must always happen, no matter how 
the parameters of the model are assigned, that v s > fbbn- 

The frequency scale v s is related to the duration of the stiff phase and it is bounded from below by 
the nucleosynthesis frequency. The slope of the spectral energy density in the high-frequency branch is 
related, ultimately, to the sound speed and it is bounded, from above, by the speed of light. These are 
the two supplementary parameters of TACDM scenario. In Fig. 4 (see plot at the left) the spectral 
energy density computed in the TACDM scenario is illustrated for two different values of wt and rx- 
For v > v B the spectral energy density acquires a blue spectrum '. 

Defining as H ~ (enA-n)^' 2 Mp the typical inflationary curvature scale and as H r the Hubble rate 



7 The spectrum is blue, in general terms, if it is increasing with frequency. The slow-roll dynamics always implies, 
within the ACDM scenario, much milder scaling violations which involve only red spectra, i.e. spectra which are very 
slowly decreasing in frequency like those of Fig. 1 (see plot at the left). 
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at the onset of the radiation epoch, u s and f ma x can be written, in the TACDM scenario, as 

al ^ 7 ^_ A 2+I ( ffiJRO \ 1/4 

4.15 x 10" 



1.177 x 10 n ST {iteVn)— ( A , [ ° , n _ 5 ) Hz, (28) 



7+1 , „ „ 1 



)lly ,.i r fr§"RQ ^ l/4 h, r = ( ±y* (Jk\ - 



zw = 1.177 x lO^S- 1 u 1V " - Hz, E = — , (29) 

V4.i5xio- 5 / ' vm p ; VM P y ' v ' 

where 7 = 7(i«t) = 3(i0t + l)/(3wt — !)• By definition E is fully determined by fixing H r . So, E and 
7(u>t) can be chosen as the two pivotal parameters of the TACDM scenario. Equivalently E and 7 can 
be traded for v 6 and for the slope of the spectral energy density during the stiff phase (which is given, 
up to logarithmic corrections) by (Qwt — 2)/[(3u?t + !)]• As stressed above, the natural upper limit for 
the spectral slope is exactly 1 which is the maximally stiff fluid compatible with causality [19]. 

The frequency v s can be much larger than zr,bn (for instance v ~ mHz in [15]) but cannot be 
smaller than Z7,t> n which constitutes a natural lower limit for v s . If u s < fbtm the pl asma would be 
stiff also throughout nucleosynthesis which is unacceptable. The observed abundances of the light 
elements (together with CMB data) also constrain the total energy density of the relic gravitons, i.e. 
the integral of Oqw^to) over the frequency. This bound is usually expressed as 8 : 

T^max / /l^O \ 

h 2 n GW ( T0 ) = h 2 O GW (^T )riln^ = 5.6xl0- 6 ° ^° )AN V (30) 

where A^ is the equivalent number of extra-relativistic species at the onset of standard big-bang 
nucleosynthesis 9 . In the standard scenario for the synthesis of light nuclei, 0.2 < AN U < 1 and, 
therefore /iof2Gw( r o) wn l be constrained accordingly. In Fig. 2 (plot at the right) the spectral energy 
density is reported as function of tt in the context of the TACDM scenario and for typical frequencies 
in the operating window of wide-band interferometers. As r^ diminishes, the amplitude of the spectral 
energy density is almost constant. The latter occurrence arises for two independent reasons. On one 
hand the most relevant constraint, in the case of growing spectral energy densities, is the one provided 
by Eq. (8) and enforced in both plots of Fig. 2. On the other hand the frequency u s depends also 
upon tt (through e, see Eqs (14) and (29)). It should be finally appreciated, from Figs. 2 and 3, that 
the pulsar timing bounds (recently revisited [32]) still imply that /iq^gw (^pulsar > r o) < 1-9 x 10~ 8 for 
a fpuisar — 10 nHz which is roughly comparable with the inverse of the observation time along which 
the pulsars timing has been monitored. Such a bound is not constraining for the TACDM model. 
The proof goes as follows. Assuming the maximal growth of the spectral energy density (i.e. that 
/igOGw(f,To) oc v) and the minimal value of v s (i.e. u > Z7,bn)> we will have that, at the frequency 



Coherently with established conventions In will denote the natural logarithm, while the logarithms to base 10 (i.e. 
common logarithms) will be denoted by log. 

9 The language of Eq. (30) may seem a bit contrived but it is a simple consequence of the historical development of 
the field. The extra-relativistic species were associated, in the past, with families of neutrinos. The nature of the bound 
on AiVj, (and hence on /lo^Gw('H))) does not change if the relativistic species are bosonic (like in the case of gravitons). 
For a discussion of the derivation of Eq. (8) see [31]. 
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scale of the pulsars, /tg^GW (^pulsar, To) — 10 13 or even 10 14 depending upon r^. But this value is 
always much smaller than the constraint stemming from pulsar timing measurements. 

In this paper it has been suggested that the ACDM parameter can be complemented by adding 
a post-inflationary phase characterized by a sound speed larger than the one of an ultra-relativistic 
plasma (i.e. 1/V3). Causality constrains the maximal barotropic index and the maximal sound speed. 
Big bang nucleosynthesis sets limits both on the maximal duration of the stiff phase and on the 
total energy density of the relic gravitons. Two new parameters will then be added to the ACDM 
paradigm which has been dubbed, throughout the paper, as tensor- ACDM (TACDM) paradigm since 
relic gravitons are copiously produced at high frequencies (i.e. larger than 0.1 nHz). The new pivot 
frequency defines the scale at which the spectral energy density of the relic gravitons starts increasing 
with a slope which is dictated by the stiff barotropic index. In the TACDM scenario, which may be 
seen as an improved version of the models proposed in [14, 15], the spectral energy density of the 
relic gravitons can even be from 6 to 7 orders of magnitude larger than in conventional inflationary 
models. Along a more technical perspective, a numerical recipe for the calculation of the spectral 
energy density has been presented. 

The advanced versions of wide-band interferometers are germane to the theme of the present 
investigation. At the moment the CMB data [1, 2], large-scale structure observations [3, 4] and 
supernovae light curves [5, 6] are used in combined analysis to put bounds on 7~t, i.e. the tensor to 
scalar ratio. Few years from now the three aforementioned cosmological data sets will still be used to 
constrain (and hopefully determine) rx while, given the foreseen sensitivities, the (terrestrial) wide- 
band interferometers will still be unable to set concurrent limits to backgrounds of relic gravitons. 
Provided the claimed sensitivities will be reached in due time, the considerations presented here give 
a concrete opportunity of using interferometers data together with the more classic cosmological data 
sets to rule out (or, more optimistically, rule in) a class of specific models. It is productive to stress 
that, in the present context, any potential upper limit from wide-band interferometers will directly 
constrain the post-inflationary thermal history. Cosmology is not tested in a laboratory: therefore 
the nature of the observations is inextricably bound to the models employed to analyze the data and 
to the potential redundancy of different data sets. It has been shown here that by complementing a 
known model with supplementary physical considerations, the three established cosmological data sets 
can also profit of a qualitatively new class of observations, such as the ones provided by wide-band 
interferometers. It is tempting to speculate that the perspective of the present paper could provoke 
a useful synergy between communities scrutinizing different branches of the graviton spectrum. The 
fruitful dialogue between the experiments sensitive to small (i.e. v v ~ aHz) and to intermediate 
frequencies (i.e. vyst — 0.1kHz) could be extended, in principle, also to conceptually different kinds of 
detectors such as microwave cavities [33] and waveguides [34]. 

It is a pleasure to acknowledge interesting discussions with E. Picasso. 
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